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Abstract 

Magnetic phase diagram of a spatially anisotropic, frustrated spin-^ Heisenberg antiferromagnet 
on a stacked square lattice is investigated using second-order spin-wave expansion. The effects 
of interlayer coupling and the spatial anisotropy on the magnetic ordering of two ordered ground 
states are explicitly studied. It is shown that with increase in next nearest neighbor frustration 
the second-order corrections play a significant role in stabilizing the magnetization. We obtain two 
ordered magnetic phases (Neel and stripe) separated by a paramagnetic disordered phase. Within 
second-order spin-wave expansion we find that the width of the disordered phase diminishes with 
increase in the interlayer coupling or with decrease in spatial anisotropy but it does not disappear. 
Our obtained phase diagram differs significantly from the phase diagram obtained using linear 
spin-wave theory. 

PACS numbers: 75.10.Jni, 75.40.Mg, 75.50.Ee, 73.43.Nq 
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I. INTRODUCTION 



Availability of new magnetic materials and the recent discovery of superconductivity at 
relatively high temperatures in the iron pnictide family of materials have spurred a flurry of 
interest in understanding the properties of frustrated magnets.-"— For the last two decades 
the properties of quantum spin-i Heisenberg antiferromagnet (HAFM) with nearest neighbor 
(NN) Ji and next nearest neighbor exchange interactions (NNN) J2 on a square lattice have 
been studied extensively by various analytical and numerical techniques.—"— It is now well 
established that at low temperatures these systems exhibit new types of magnetic order and 
novel quantum phases.— For J2 = the ground state is antiferromagnetically ordered at 
low temperatures. Addition of NNN interactions induces a strong frustration and break 
the antiferromagnetic (AF) order. The competition between NN and NNN interactions for 
the square lattice is characterized by the frustration parameter t] = J2/ Ji- A disordered 
paramagnetic phase probably columnar dimer exists between rjic ~ 0.38 and ri2c ~ 0.60. 
For 7] < rjic the square lattice is AF-ordered whereas for 77 > t]2c a degenerate collinear 
antiferromagnetic (CAF) stripe phase emerges. Experimentally by applying high pressure 
the ground state phase diagram of these frustrated spin systems can be explored from low 
?7 = J2/J1 to high 7]. For example in Li2VOSi04 X-ray diffraction measurements show that 
the value of 1] decreases by about 40% with increase in pressure from zero to 7.6 GPa.— 
Moreover, nuclear magnetic resonance, magnetization, specific heat, and muon spin rotation 
measurements on the compounds Li2VOSi04, Li2VOGe04, VOM0O4, and BaCdVO(P04)2 
show significant coupling between NN and NNN neighbors.-"- In addition these experiments 
on Li2VOSi04 have shown that it undergoes a phase transition at a low temperature (2.8 K) 
to collinear AF order with magnetic moments lying in the a — b plane with J2 + Ji ~ 8.2(1) 
K and J2/J1 ~ 

A generalization of the frustrated Ji — J2 model is the Ji — J[ — J2 model where C = J'l/ Ji 
is the directional anisotropy parameter.— "^"^ A possible candidate of this model may be the 
compound (NO)Cu(N03)3.— Extensive band structure calculations^^ for the vanadium phos- 
phate compounds Pb2VO(P04)2, SrZnVO(P04)2, BaZnVO(P04)2, BaCdVO(P04)2 have 
shown four different exchange couplings: Ji and J[ between the NN and J2 and between 
NNN. For example C, ~ 0.7 and J2/ J2 ~ 0.4 were obtained for SrZnVO(P04)2. Recently us- 
ing second-order spin-wave expansion the effects of directional anisotropy on the spin-wave 
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energy dispersion, renormalized spin-wave velocities, and magnetizations for the two ordered 
phases have been studied in detail.— It has been found that the spatial anisotropy reduces 
the width of the disordered phase. 

Although much efforts have been made to understand the properties of two dimensional 
(2D) frustrated magnets research on three dimensional (3D) systems has been limited.—"— 
Earlier work on HAFM on the pyrochlore lattice^ and on the stacked kagome lattice^"— 
showed existence of a magnetically disordered phase. On the contrary, spin-wave calculations 
for the 3D Ji — J2 model on the body-centered cubic (BCC) lattice and for the simple cubic 
lattice (SC) show no signs of an intermediate disordered paramagnetic phase.— 1^1^ Instead 
for the BCC and the SC lattice a direct first order phase transition occurs from a two- 
sublattice Neel ordered AF phase for small J2 to a collinear antiferromagnetic ordered state 
for large 72.—*^ For Li2VOSi04, a layered material that can be described by a square lattice 
Ji — J2 model with large J2 the interlayer coupling J_l/Ji ~ 0.07 is not negligible. Due to 
a finite interlayer magnetic coupling these experimental systems are quasi 2D. 

Another example is the recently discovered iron pnictide superconductors.— The parent 
phases of these materials have been found to be metallic but with AF order and magnetic 
excitations have shown to play an important role in the superconducting state.— "—'^"— 
Although magnetism in these materials are still debated neutron scattering spectra for the 
pnictides show sharp spin-waves. ^'^ These studies have also revealed that the parent com- 
pounds exhibit a columnar antiferromagnetic ordering with a staggered magnetic moment 
of (0.3 - 0.4)/iB in LaOFeAs and (0.8 - l.Ol)/^^ in Sr(Ba,Ca)Fe2As2.— >^ Moreover, at low 
temperatures there is orthorhombic distortion and the exchange constants have been found 
to be substantially anisotropics^'^- Motivated by the observation of spatially anisotropic ex- 
change constants in these materials the spin-wave spectra and the low-temperature phases 
of the model can be studied by the spatially anisotropic Ji — .][ — J2 Heisenberg model 
on a square lattice with NN exchanges Ji along the x axis, J[ along the y axis, and NNN 
interactions J2 along the diagonals in the xy plane. Recent experiments on iron-based super- 
conductors such as undoped iron oxypnictides reveal that the electronic couplings are more 
three dimensional than in the cuprate superconductors.— 1^>^ With decrease in temperature 
most undoped iron-pnictide superconductors show a structural transition from a tetragonal 
paramagnetic phase to a orthorhombic phase. In the 122 materials a three dimensional 
long-range antiferromagnetic order develops simultaneously. 
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These quasi 2D frustrated layered systems with anisotropic magnetic exchange couphngs 
and with non-neghgible interlayer couphngs J±_ serve as a motivation to investigate the Ji — 
J[ — J2 — J± modeL This model on a square lattice with isotropic NN interactions ( Ji = J[) 
has been studied numerically by coupled-cluster and the rotation-invariant Green's function 
methods." These calculations show that the quantum paramagnetic phase disappears for 
a critical value of interlayer coupling JI ^ 0.23Ji. For J± < JI a second-order phase 
transition between the quantum Neel to a quantum paramagnetic phase and then a first- 
order transition from the quantum paramagnetic phase to the stripe phase occur. For J± > 
JI there is a direct first-order transition between the Neel to the stripe phase. Existence of 
this critical point has also been reported by other authors where they have used effective 
field theory in a finite cluster to obtain JI ^ 0.67Ji.— In the context of iron pnictides 
this model has been studied using linear spin-wave theory (LSWT) to obtain the spin-wave 
energy dispersion and the sublattice magnetization.— However, to our knowledge this model 
has not been studied using higher order spin-wave expansion. We will find that second- 
order corrections due to quantum fluctuations increase substantially as the classical phase 
transition point is approached. As a result the magnetic phase diagram obtained from 
second-order spin-wave expansion differs significantly from the phase diagram obtained by 
LSWT. 

In this work we investigate the magnetic phase diagram of the Ji — J{ — J2 — J± Heisenberg 
AF on a stacked square lattice. We use spin-wave expansion based on Holstein-Primakoff 
transformation up to second-order in 1/S to numerically calculate the sublattice magne- 
tization for each of the two ordered magnetic phases. The paper is organized as follows. 
Section [TTl provides an introduction to the Hamiltonian for the Heisenberg spin-i AF on a 
spatially anisotropic stacked square lattice. The classical ground state configurations of the 
model and the two phases are then briefly discussed. In the next two sections Sec. Ill Al 
and Sec. Ill Bl the spin Hamiltonian is mapped to the Hamiltonian of interacting spin- wave 
excitations and spin-wave expansion for sublattice magnetizations are presented for the two 
phases. Magnetizations for the two phases are numerically calculated with different values 
of interlayer coupling and spatial anisotropy and the results are plotted and discussed in 
Section IIIII Finally we summarize our results in Section IIVI 



4 



II. THEORY 



We consider a spatially anisotropic, frustrated spin-| HAFM on a x A^^ x cubic 
lattice with four types of exchange interactions between spins: Ji along the x (row) direction, 
J[ along the y (column) direction, J2 along the diagonals in the xy plane, and J± is the 
interlayer coupling. We assume all interactions to be AF and positive i.e. Ji, J(, J2, J± > 0. 
This Ji — J[ — J2 ~ J± spin system is described by the Heisenberg Hamiltonian 



H = - 



where ^ labels the layers, i runs over all lattice sites and i + S^ {^x = =tl) and i + 5y {6y = ±1) 
are the nearest neighbors to the i-th site along the row and the column direction. The third 
term represents the interaction between the next-nearest neighbors, which are along the 
diagonals in the xy plane and the last term is for the NN coupling between the layers. This 
model is different from the fully frustrated simple cubic lattice which has the additional 
NNN interactions in the xz and yz planes.— 

At zero temperature this model exhibits three types of classical ground state (GS) con- 
figurations: the Neel or the (tt, tt, it) state and the two stripe states which are the columnar 
stripe (vr, 0, vr) and the row stripe (0, vr, vr). The classical ground state energies of these states 
are 

Efr/N = -^JiS'z [1 + C - 2r7 - 5] , 

Ef^/N = -\j,S^z [1 - C + 2r/ + 5] , (2) 

Ef^VN = -^JiS'z [-1 + C + 2r/ + 5] , 

where ( = J[/ Ji measures the directional anisotropy, = J2/J1 is the magnetic frustration 
between the NN (row direction) and NNN spins, and 6 = Jj_/ Ji is the interlayer coupling 
parameter. 2; = 2 is the number of NN sites along the row (column) direction. For t] < (/2 
the classical GS is the AF Neel state and for rj > 1/2 (with C = 1) the GS is doubly 
degenerate. Otherwise for C < 1 the GS is the columnar AF (vr, 0, vr) state. The classical 
first-order phase transition between the AF and CAF state occurs at the critical value 
j^ciass _ ^^2^ which is independent of 6. 

The motivation of this paper is to investigate the role of interlayer coupling S and spatial 
anisotropy ( to the quantum phases of this model. How does the quantum fluctuations 



due to 5 and ( affect the disordered paramagnetic phase? Is there a critical point for S 
(or C) above (or below) which the intermediate quantum paramagnetic GS does not exist 
and we have a direct first-order phase transition from the AF to the CAF ordered phase? 
For our study we follow the standard procedure by first expressing the fluctuations around 
the "classical" ground state in terms of the boson operators using the Holstein-Primakoff 
transformation. The quadratic term in the boson operators corresponds to the LSWT, 
whereas the higher-order terms represent spin-wave (magnon) interactions. We keep terms 
up to second order in 1/S. The staggered magnetization per spin to the leading order in 
l/S*^ for the AF and CAF phases are then obtained from the renormalized magnon Green's 
functions and self-energies. We will follow the theoretical framework described in detail in 



Ref. 



23l . However for completeness we provide the necessary equations that are required for 



numerical computations and to follow the present work. 



A. (vr, vr, n) AF Neel Phase 

For the AF ordered phase NN interactions are between A and B sublattices and NNN 
interactions are between A-A and B-B sublattices. The Hamiltonian in Eq. [1] takes the form: 

H = Ji ^ ■ S^5^ £ + J[ ^ S^^ ■ S^5^^^ + -J2 ^ |^S^£ ■ S^5^+5^ £ + ■ S|^5^+5^_£ 

i,l i/ i,e 

+ J±^Sf'e ■ S|^^+i. (3) 

This Hamiltonian is mapped into an equivalent Hamiltonian of interacting bosons by trans- 
forming the spin operators to bosonic creation and annihilation operators a\ a for "up" and 
b'^ , b for "down" sublattices using the Holstein-Primakoff transformations keeping only terms 
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up to the order of l/S*^ 



qA+ 



qA- ^ 



an, 



2Sa 



1 alfCiie 1 dl^diialidii' 
~ 2 (25) ~ 8 (25)2 . 

2 (25) ~ 8 (25)2 



1 - - 



5, 



5 — al^aa, 



qB+ 



2Sb], 



1 - - 



1 bl.bji 1 h]fbjth]fbjt 



S 



1 - - 



2 (25) 8 (25)2 
2 (25) 8 (25)2 



In powers of 1/5 the Hamiltonian is now written as 



H 



1 



1 - 



2r] 



c 



+ H0 + HI + H2 + .... 



(4) 



(5) 



The first term corresponds to the classical energy of the AF ground state (Eq. [2]). Then the 
real space Hamiltonian is transformed to the k-space Hamiltonian. Momentum k is defined 
in the first Brillouin zone (BZ): — vr < < 71, — vr < ky < tt, — tt < k^ < ti (with unit 
lattice spacings). Next we diagonalize the quadratic part Ho by transforming the operators 
Ok and 6k to magnon operators and /3k using the Bogoliubov (BG) transformations 



al = l^al + mk/3-k, &-k = "^k"! + ^k/3-k, 



(6) 



where the coefficients /k and mk are defined as 



with 



Ck 



2ek 



1/2 



mk 



-sgn(7i, 



Ck 



2ek 



1/2 



— Xk/l( 



(7) 



Ck 

7k 

7lk 
72k 

Kk 



(1 - 7^)^/^ 

7ik/«;k, 



cos(A;a;) + (cos{ky) + 5cos(A;2)]/ (1 + C), 

5 



cos(A;2:) cos(A;j^), 
2r] 



1 - 



1 + C 



1 - 72k) + 



(8) 
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7k is negative in certain parts of the first BZ - so it is essential to keep track of the sign 
of 7k through the function sgn(7k). After these transformations, the quadratic part of the 
Hamiltonian becomes 

Ho = J,Sz{l + C) I] /^k (ek - 1) + JiSz{l + C)Y1 ^k^k («i.ak + /^^/^k) • (9) 

k k 

The first term is the zero-point energy and the second term represents the excitation energy 
of the magnons within LSWT. 

The part Hi corresponds to 1/5* correction to the Hamiltonian. We follow the same pro- 
cedure as described above. The resulting expression after transforming the bosonic operators 
to the magnon operators is 

JiSz{l + C) 

k 



Hi 



2S 



2SN 



1234 



+ h.c. 



(10) 



In the above equation three-dimensional momenta ki,k2,k3,k4 are abbreviated as 1, 2, 3, 
and 4. The first term in Eq. [10] is obtained by setting the products of four boson operators 
into normal ordered forms with respect to the magnon operators, where Ak and i?k are 



with 



Ak 
5k 



Av 



Bi 7ik 

/tkCk 



«^k - 7ik 



+ A2 — 

Ck 



1 - 72k 



1 - 72k 



(11) 

(12) 



A2 



N 



p L ft^p 



^P-I 



27] 



l + CJ N 



- 5^ ^ 1 1 - ep - 72p 



f 27] \ 2 ^ 1 r 7i'pl 



(13) 
(14) 

(15) 



The second term in Eq. [10] represents scattering between spin-waves where the three- 
dimensional delta function 5g(1 + 2 — 3 — 4) ensures that momentum is conserved within a 
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reciprocal lattice vector G. Explicit forms of the vertex factors ^1*234 "^ are given in Ref. |23 



Here we provide two of the vertex factors that are needed to calculate the magnetization. 
They are 

^1234 = -71(2 - 4)3:4 - 7l(l - 4:)xiX2X4 - 7i(2 - 3)X3 - 7i(l - 3)xiX2Xs 

+ ^ 7i(2) + 7i(l)a;ia;2 + 7i (3)^2X3 + 7i (4)0:2X4 

+ 7i(2 - 3 - 4)0:3X4 + 7i(l - 3 - 4)xiX2X3X4 + 7i(3 - 2 - l)xiX3 + 71 (4 - 2 - l)xiX4 

27] 



1 + C 



/: 



1234 



X2 + Sgn(7G)XiX3X4 



(16) 



(6) 



1234 



7i(2 - 4)x2X3 + 7i(2 - 3)x2X4 + 71(1 - 3)xiX4 + 71(1 - 4)xiX3 



- - 7i(2)x2X3X4 + 7i(3)x4 + 7i(2 - 3 - 4)x2 + 71 (3 - 2 - l)xiX2X4 

+ 7i(l)xiX3X4 + 7i(4)x3 + 7i(l - 3 - 4)xi + 71 (4 - 2 - l)xiX2X3 
2r] 



1 + C 



1234 



X3X4 + Sgn(7G)xiX2 



(17) 



with 



1 r 



1234 



72(l-3) + 72(l-4) + 72(2-3)+72(2-4)-72(l)-72(2)-72(3)-72(4) . (18) 



The second order term, H2 is composed of six boson operators. After transformation to 
magnon operators ak, /3k the Hamiltonian in normal ordered form reduces to 
JiSzil + 



Ho 



{2sy 



J2 (alftk + I^IM + C2k (ttl/3lk + /3-k«k) + • • • 



(19) 



The dotted terms contribute to higher than second-order corrections and are thus omitted 
in our calculations. The coefficients Cik and C2k are 



Cik = ^1(^7) X]^i^2 - 671(2 - 1 - k)xkXiX2 + 7i(2)xiX2 + 7i(2)xkXiX2 



12 

+ 27i(k)xkXi + 7i(l)XkXi + 7i(2)x2 
2^Kn) 



2r] 
4VlTC 



ll{l+xl)C^, 



(20) 



C2k = ^1(^7) XI ^1^2 371(2 - 1 - k)xiX2 + 371(2 - 1 - k)XkXiX2 - 27i(l)xkXix: 



12 



27i(2)xkX2 - 7i(k)x2 - 7i(k)xkX2 



27] 

2VlTC 



^k'^kC'k, 



(21) 



with 



Ck 



2 \2 



N 



f272(k) + 72(1) + 72(2) - 472(k +1-2) 



12 



xl 



+ 



72(2)-72(l + 2-k)](l + x?x^)}. 



(22) 



After defining the renormalized Green's function as in Ref. |23| the first and second order 
self-energies are written as 



CO 



S^^j(k,a;) = Ak, 
S(^i(k,a;) = 5k, 



(23) 
(24) 



pq 



X 



|y(4) 1 2 

I k,p,q,[k+p-q] I 



|T/{6) 12 
I k,p,q,[k+p-q]l 



.U — Ep — Eci — -Ek+p-q + U + Ep + + -Ek+p-q 



iS. 



(25) 



sS(k,a;) 



Si'j,f-k, -w) 



2k 



X v; 



(4) 



X. 



(6) 



pq 

2{Ep + i^q + -Ek+p-q) 



(26) 



k,p,q,[k+p-q] k,p,q,[k+p-q]^2 _ + E^ + ^k+p-q)^ ' 

where [k + p — q] is mapped to (k + p — q) in the first BZ by the reciprocal vector G. 

The magnetization M defined as the average of the spin operator Sz on a given sublattice 
(say A) is expressed as 



M = S-{ala,) = S-AS + ^+ 



{2S) {2Sy' 



(27) 



where 



Ml 
M2 



k 

2 

iV 



(28) 
(29) 



'-'*ES(k.-B.) 



^^V^ 2/2/2/2/2 



^k 



k,p,q,[k+p-q]l 



pq 



2/kmkSgn(7G)Vj^^'^q^,^^p_qjV;^pqj;^^p_qj 

-^k ~ (-^P + -^q + -£^k+p-q)^ 



(Ek + -£^p + -E'q + -Ek+p-q)' 
K6) 



+ 



(30) 



The zeroth-order term AS* corresponds to the reduction of magnetization within LSWT, Mi 
term corresponds to the first-order 1/S correction, and M2 is the second-order correction. 
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B. (7r,0,7r) CAF Phase 



The Hamiltonian describing the CAF phase is 



1 



£+1- 



(31) 



The procedure is same as the AF phase. For this phase the structure factors 7^^,^, along 
with other quantities required for the calculations are defined as 



7^^ = 


[cos(/c^)(l 


72k = 


cos(A;y), 


7k = 


7ik/<> 


/ 

Kjj. — 


1- ^ 

1 + 27] 


4 = 





7^ + 



1 + 27]' 



(32) 



The coefficients that appear in the Hamiltonian Hi are 



with 



1 



/2 



7ik 



CkL 



72k 



'^k^k 



72k 



/2 



^2 - (tt2;^)|5Zz[^ 



27] 



72r 



E;rf72p 



p P 



1% 



(33) 
(34) 



(35) 
(36) 
(37) 



ifo) Hi, and can be expressed in the same forms as in Eqs. [9], [101 ^"^^ [19] with the new 
coefficients AJ^, i?^, CJj^, C2J5. and with the replacement o 27/. The expressions for the two 
vertex factors V'^^\V''^^'^ and the coefficients C^^jCgk are similar to the AF phase (details 



can be found in Ref. 



23|). 
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III. MAGNETIZATION AND THE PHASE DIAGRAM 



We obtain the sublattice magnetization M for the two ordered phases with different values 
of (, 1], and 6 from Eq. [27]by numerically evaluating Eqs. I2BH5U] (using similar expressions for 
the CAF phase). Especially to obtain the second order correction term M2 we sum up the 
values of N^/S points of k in the 1/8-th part of the first BZ and points of p and points 
of q in the full BZ. To check the convergence of our results we do the calculations for the 
AF-phase with Nl = 8, 10, and 12 sites for C = 1 and S = 0.1. The convergence is very good 
as shown in Fig. [TJ Maf becomes zero at the critical point rji^ ~ 0.460. Hereafter, we use 
Nl = 12 lattice sites for all of our numerical computations. Evaluation of the magnetization 
requires summing contributions from over 645 million points in the first BZ for each (, 1], 
and 6. 

0.4 
0.3 
^< 0.2 
0.1 


0.1 0.2 0.3 0.4 0.5 0.6 

Tl 

FIG. 1: (Color online) Convergence test of magnetization calculation for the AF phase (C = 1;<J = 
0.1) with Nl = 8 (dashed line), 10 (solid line), and 12 (open circles) lattice sites . The convergence 
is excellent with Nl = 12 sites. 

Figure [2] shows the sublattice magnetization M^^ for the AF phase with C = 1 ^ind 
f] = for different values of interlayer coupling 6. In the inset we plot the spin-deviation 
A = 0.5 — AI^^ with 6. A is a measure of quantum fluctuations from the classical value of 
0.5. We find that with increase in 6 the fluctuations decrease, thus M^^ (with second-order 
corrections) increases from 0.308 for 5 = to 0.423 for 6 = 1.0. This result is expected 
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1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

C=1, 6=0.1 1 


1 1 1 1 1 1 1 


Nl=8 ; 




— Nl=10 




oNl=12 






) 





\ 




as with the increase in interlayer couphng the system undergoes a dimensional transition 
from 2D to 3D. The dotted hne is the result from LSWT. Result for Mq^ obtained from 
LSWT captures the essential physics both qualitatively and quantitatively as second-order 
corrections are small for the unfrustrated (r/ = 0) case. 




FIG. 2: (Color online) Magnetization M^^ for the AF phase plotted for the unfrustrated case 
(?7 = 0) with interlayer coupling 5 for C, = 1.0. Dotted line is the result from LSWT and the 
solid line is with l/S"^ corrections. (with second-order corrections) increases from 0.308 for 

5 = to 0.423 for 5 = 1.0. Inset shows the spin-deviation A = 0.5 — M^^ . With increase in 6 the 
fluctuations decrease, thus Mq^ increases. This result is expected as with increase in the interlayer 
coupling the system makes a transition from two to three dimensions. 

Figures |3] and H] show the magnetization with increase in the frustration parameter 
for two different values of the spatial anisotropy parameter = 0.9 and 0.6. Similar to 



the 2D case (see Ref. 



23 



for details) our spin-wave expansion for the CAP phase becomes 
unreliable as ( gets close to 1. Thus we have not chosen ( = I for our plots. For each ( 
and 6 three different curves are plotted: the long-dashed lines represent LSWT results, the 
dotted lines include the first-order 1/S corrections, and the solid lines represent corrections 
up to second-order to the LSWT results. As t] approaches the classical transition point 
^ciass _ g_^g from both sides of the two ordered phases the dotted curves diverge. However, 
l/S"^ corrections (M2) significantly increase to stabilize the divergence. We find that the 
magnetizations with 1/5*^ corrections decrease steadily and then sharply drops to zero for 
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both the phases as 77 — )■ rif'^^^. As an example, with (^ = 0.9, 5 = 0.1 magnetization for the 
AF phase begins from ^ 0.377 at r/ = 0, then decreases till 1] ^ 0.41, and sharply becomes 
zero at the critical point rji^ ~ 0.427. M2 corrections start from a small positive number at 
rj = and then switches sign at rj ^ 0.32. M for the CAF phase shows the same feature, 
where we find M ^ 0.393 at 77 = 1 and r]2c ~ 0.464. However, in this case M2 corrections are 
always negative. These results are qualitatively similar to the results for the 2D spatially 
anisotropic, frustrated HAFM on a square lattice {6 = 0).— 

On the other hand, LSWT calculations for ( = 0.9 and 6 = 0.1 show that the magne- 
tization becomes zero for the AF phase at r/ic ~ 0.45 whereas it does not go to zero from 
the CAF phase. LSWT is not applicable at the classical transition point rj = 0.45. Extrap- 
olation of the CAF phase LSWT results show a direct transition from the AF to the CAF 
phase. 

For ( = 0.9 and S = 1.0, M never becomes zero for both the AF and CAF ordered phases 
(within LSWT) indicating that there is a direct first-order phase transition from the AF to 
the CAF phase. But with second-order corrections we find that M becomes zero for both 
the phases at the critical points rjic ~ 0.434 and ri2c ~ 0.469. With increase in 6 the width 
of the disordered PM phase diminishes but it survives even for large 6. 

In Fig. Hlwe show the results for ( = 0.6. LSWT calculations show a first-order direct 
phase transition from AF to the CAF phase for both 6 = 0.1 and 1.0. Similar to the ( = 0.9 
case we find that M vanishes (with second-order corrections) for both the phases at the 
critical points rjic ^ 0.288, r^ac ^ 0.311 for 6 = 0.1 and r/ic ^ 0.290, r^ac ^ 0.310 for 6 = 1.0. 
The two phases are separated by a narrow disordered paramagnetic region. We repeat the 
calculations for ( = 0.4 with different values of S and obtain similar features (the results are 
not shown). 

Another feature we observe from our data is that with decrease in spatial anisotropy ( the 
width of the disordered region (772c ~ Vic) diminishes but it never disappears. This is shown 
in Fig. [5] for 6 = 0.1. The solid lines represent the critical points rjic and 772c for the AF and 
CAF phases. The dashed line is the classical first-order phase transition line 77^^^'^'^ = C/2 
(independent of 6) between the two phases. For the AF phase rjic ~ 0.460 for C = 1- For 
the CAF phase we obtain the value of the critical phase transition point ri2c ~ 0.53 at = 1 
by extrapolation as our spin- wave expansion is unreliable for ( near 1. These results can be 
compared to the results for the spin-1 spatially anisotropic, frustrated HAFM on a square 
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FIG. 3: (Color online) Sublattice magnetization M for the AF and CAF ordered phases for spatial 
anisotropy C, = 0.9 with two different values of interlayer couplings 5 = 0.1 and 1.0. For each 5 
results from LSWT (long-dashed lines), with first-order 1/S corrections (dotted lines), and with 
second-order l/S"^ corrections (solid lines) are plotted. With increase in rj dotted lines diverge. 
Second-order corrections become significant with increase in r] and stabilize the anomalous diver- 
gence of the magnetization. LSWT calculations for 5 = 0.1 show that the magnetization becomes 
zero for the AF phase at ??ic ~ 0.45 whereas it does not go to zero from the CAF phase. LSWT is 
not applicable at the classical transition point rf^^^ = 0.45. Extrapolation of the CAF phase LSWT 
results show a direct transition from the AF to the CAF phase. For 6 = 1.0, M never becomes zero 
for both the phases (within LSWT) indicating that there is a direct first-order phase transition 
from the AF to the CAF phase. However, with second-order corrections M becomes zero for both 
the phases at the critical points r/ic ~ 0.427, r/2c ~ 0.464 for 5 = 0.1 and rjic ~ 0.434, ?72c ~ 0.469 
for 5 = 1.0. 



lattice {6 = 0) where we have found t]ic ~ 0.41 and 772c ~ 0.58 for ( = 1.^^. All other features 
are similar to the phase diagram of the anisotropic, frustrated 2D model. 

Phase diagram for the Ji — J[ — J2 — J± model for C = 0.9 is displayed in Fig. |6l Within 
our spin-wave expansion we do not find any critical value of Jj above which there is no 



disordered region. This is in contrary to the findings in Refs. 



72 



and 



73 



Instead from our 



calculations we find that the transition between the two phases is always separated by the 
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C=o.6, 5=1.0 ; 

-- LSWT 

with 1/S 

— with 1/s' 

0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

FIG. 4: (Color online) Magnetization M for the AF and CAF ordered phases shown for spatial 
anisotropy C, = 0.6 with interlayer couplings 5 = 0.1 and 1.0. Similar to Fig. [3] M with 1/ S"^ 
corrections decrease steadily and then sharply drops to zero for both the phases. LSWT calculations 
for both 5 =0.1 and 1.0 show that there is a direct first-order phase transition from the AF to the 
CAF phase. However, with second-order corrections we find that M vanishes for both the phases 
at the critical points r/ic ~ 0.288, r]2c ~ 0.311 for 5 = 0.1 and 7?ic ~ 0.290, 772c ~ 0.310 for 5 = 1.0. 

magnetically disordered phase. For the AF phase rjic increases from 0.386 for 5 = to 
0.434 for 5 = 1.0. On the other hand for the CAF phase ri2c begins from ^ 0.484 and then 
slowly decreases to ~ 0.462 for 6 = 0.3. Then from 6 ~ 0.4 the phase boundary for the CAF 
phase shows a slight upward rise. In this strong interlayer coupling limit where J± becomes 
comparable to Ji and J[ for the CAF phase (J2 > Ji) our model becomes unrealistic as 
we have excluded the NNN interactions in the xz and yz planes. Our present model differs 
from the spin-| Heisenberg AF on a fully frustrated simple cubic lattice where we have the 
additional J2 interactions in the xz and yz planes. In that case it has been shown that there 
is a direct first-order phase transition from the AF to the CAF phase.— 
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FIG. 5: (Color online) Phase diagram for the Ji — J'l — J2 — J± model with 6 = 0.1. Solid lines 
represent the critical phase transition points ??ic and r]2c for the AF and CAF phases. Dashed line 
shows the classical phase transition line rjf'^'^'^ = (^/2. For the AF phase rjic ~ 0.460 at C = 1. 
Spin- wave expansion for the CAF phase becomes unreliable for ^ near 1. We thus extrapolate our 
data to obtain r/2c ~ 0.53 for C = 1- The width of the disordered region (r/2c — Vic) increases with 
the anisotropy parameter ^. 

IV. CONCLUSIONS 

In this work using second-order spin-wave expansion we have studied the effects of in- 
terlayer coupling and directional anisotropy on the magnetic phase diagram of a frustrated 
spin- 1 Heisenberg antiferromagnet on a stacked square lattice. Linear spin- wave theory 
calculations for this model show that for small interlayer coupling there are two magnetic 
ordered phases, AF and CAF which are separated by a disordered paramagnetic state. How- 
ever when the interlayer coupling exceeds a critical value the disordered paramagnetic phase 
disappears and then there is a direct first-order phase transition from the AF to the CAF 
phase. Recent numerical calculations using coupled cluster and rotation-invariant Green's 
function methods support this picture.—*^ With our second-order spin-wave expansion we 
have found that with increase in next nearest neighbor frustration, l/S"^ corrections play 
a significant role in stabilizing the magnetization as the classical phase transition point is 
approached. As expected from linear spin-wave theory we have found that there are two 
ordered magnetic phases (Neel and stripe) which are separated by a paramagnetic disor- 
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FIG. 6: (Color online) Phase diagram for the Ji — J'l — J2 — J± model with = 0.9. The dashed 
line shows the classical phase transition line r]^'^^^ = 0.45. Phase boundaries of the two ordered 
phases AF and CAF diminishes but never disappears with increase in 6. These two ordered phases 
are always separated by the magnetically disordered phase. For the AF phase rjic increases from 
0.386 for 5 = to w 0.434 for 6 = 1.0. Whereas for the CAF phase rj2c begins from 0.484 
and then slowly decreases to ~ 0.462 for 6 = 0.3. But from 5 ~ 0.4 the phase boundary for the 
CAF phase shows a slight upward rise. By excluding the J2 interactions in the xz and yz planes 
our model becomes unrealistic in the strong interlayer coupling limit where J± is comparable to 
Ji and J[ (especially in the CAF phase). 

dered phase. But the values of the critical phase transition points for the two phases differ 
from the LSWT predictions. Our calculations show that the width of the disordered region 
diminishes with decrease in the directional anisotropy. These features are similar to the mag- 
netic phase diagram of a two-dimensional frustrated Heisenberg spin-| antiferromagnet.— 
However, with increase in the interlayer coupling we have found that the parameter region 
of this disordered phase does not disappear. Our obtained phase diagram is significantly 
different from the phase diagram obtained using linear spin-wave theory which predicts a 
direct first order phase transition from the AF to the CAF phase beyond a critical value of 
interlayer coupling. In summary with our present approach based on second-order spin wave 
expansion we do not find existence of any critical interlayer coupling (or spatial anisotropy) 
beyond (or below) which there is a direct transition from one phase to the other ordered 
phase. 
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